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Abstract 

Wc aim to assess the effectiveness of feedback controlled resonant drift pacing as a method 
for low energy defibrillation. 

Antitachycardia pacing is the only low energy defibrillation approach to have gained 
clinical significance, but it is still suboptimal. Low energy defibrillation would avoid adverse 
side effects associated with high voltage shocks and allow the application of ICD therapy 
where it is not tolerated today. 

We present results of computer simulations of a bidomain model of cardiac tissue with 
human atrial ionic kinetics. Re-entry was initiated and low energy shocks were applied with 
the same period as the re-entry, using feedback to maintain resonance. We demonstrate 
that such stimulation can move the core of re-entrant patterns, in the direction depending 
on location of electrodes and a time delay in the feedback. Termination of re-entry is 
achieved with shock strength one order of magnitude weaker than in conventional single- 
shock defibrillation. 

We conclude that resonant drift pacing can terminate re-entry at a fraction of the shock 
strength currently used for defibrillation and can potentially work where antitachycardia 
pacing fails, due to the feedback mechanisms. Success depends on a number of details which 
these numerical simulations have uncovered. 

Keywords Re-entry; Bidomain model; Resonant drift; ICD; Defibrillation; Antitachycar- 
dia pacing; Feedback. 
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Introduction 

Several clinical trials established that the timely application of an electric shock, particularly 
with ICD, is the only reliable therapy to prevent sudden cardiac death |l|). However, the strong 
shocks required are reported to have serious adverse effects, most prominently via electropora- 
tion, alterations of the action potential waveform and duration depolarization of the resting 
potential (j3|), increased pacing thresholds a, loss of excitablhty fl, and transient ectopic 
afterdepolarizations which may initiate postshock arrhythmias (0). Other studies reported me- 
chanical aftereffects such as mechanical dysfunction (stunning), increases in contractility (0) 
and hemodynamically mediated symptoms (0). Psychological effects on patients play a non- 
negligible role ([^), and despite long-term survival benefits, patients suffering from arrhythmias 



which are not immediately life-threatening do not tolerate ICD therapy (lid ). Diphasic (jlll ). 
multiphasic (fl^ ). and truncated exponential (0) shock waveforms defibrillate at a lower thresh- 
old, but still too high for painless defibrillation. 

Several approaches to minimize defibrillation energy by employing smarter protocols are 
under examination, but so far only antitachycardia pacing has gained clinical significance. An- 
titachycardia pacing is a series of weak shocks applied at a frequency higher than the intrinsic 
frequency of the arrhythmia. This therapy has conventionally been applied to slower, presum- 
ably hemodynamically tolerated, ventricular tachycardias. Fast ventricular tachycardias (188 
to 250 beats/min) typically receive high amplitude shock therapy, even though antitachycardia 
pacing may work (Il4l ). Although the mechanisms responsible for antitachycardia pacing failure 
are not fully understood, an inherent weakness is evident: a fixed pacing frequency is likely to 
be suboptimal depending on the arrhythmia. No arrhythmia-specific input is used to form or 
adjust the antitachycardia pacing sequence. 

In this study we assume that arrhythmias are sustained by re-entry and consider a method 
which employs feedback-driven pacing to control and eliminate the re-entry cores, by moving 
them until they hit inexcitable obstacles or each other, and annihilate. The method relies on 



a phenomenon of "resonant drift" (jl5l . Il6l ): the drift of re-entrant waves when periodic, low 
energy, shocks are applied in resonance with the period of the re-entry. A feedback algorithm 
(firl ) is used to maintain the resonance. Resonant drift and its feedback control, herein referred 
to as resonant drift pacing, have only been studied experimentally in the Delousov-Zhabotinsky 
reaction medium and in simulations of simplified models of cardiac tissue {e.g. (18, H, 2ol)). 



The goal of the present simulation study is to investigate the effectiveness of using resonant 
drift pacing for low- voltage defibrillation. Previous studies of resonant drift were in models which 
were very different from modern models of cardiac tissue in many important respects, and their 
relevance for low-voltage defibrillation is debatable. Here we use an anisotropic bidomain model 
of cardiac tissue with microscopic heterogeneities and realistic cellular kinetics of human atria. 
Shocks are applied by injecting current into, and withdrawing from, the extracellular space. 
Such description was used in simulations before, but only to study high-voltage single-shock 
defibrillation ^ 

Our results show that in this model setting, resonant drift pacing can be used to move 
the core of re-entrant activation patterns. Termination can be achieved with high probability, 
and within the time deemed acceptable by clinicians for antitachycardia pacing to work, at a 
fraction of the conventional single shock defibrillation strength, by moving the cores until they 
hit an anatomical boundary, or annihilate with each other. We show that in a realistically 
anisotropic model, direction of movement depends on electrode location and time delays of 



the shock application, which is in agreement with previous studies e.g. (|l5l. Il9l). Knowing the 



electrode location and the anatomy of the heart, the best delay could be estimated to move the 
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core in the direction of a suitable anatomical structure or boundary (i.e. an inexcitable piece of 
tissue) which is likely to terminate the re-entrant wave. 

Methods 

Governing equations 

The bidomain model of cardiac tissue is most widely used to study defibrillation-related phe- 
nomena 113). The system can be written as 

dV I I 

Cm-Q^ = -hon + ^y ■{cTiVV) + -V ■{cTiV(t>e), (1) 

V • ((cTi + <Te) V0e) = -V-{a,VV)+h, (2) 

where (pi and (f)e are the intracellular and extracellular potential distributions, (i is the average 
cell's surface to volume ratio, Cm is the membrane capacitance per unit area, y = c/jj — 0e is the 
transmembrane potential, (Tj and (Te are the intracellular and extracellular conductivity tensors 
respectively, Jg is the extracellular current, and lion is the ionic current density through the 
membrane. 

In absence of extracellular current Jg, if variations of (pe are negligible compared to or 
if the anisotropy ratios in both domains are the same, then ([T|) and Q can be replaced by a 
simpler monodomain equation, 

dV 1 

= -hon + ^y ■{cTmVV), (3) 

where am is the monodomain conductivity tensor. 

We used a bidomain description of cardiac tissue, with the exception of monodomain de- 
scription later than 10 ms after shocks when determining single-shock defibrillation thresholds. 
A recent comparison of the monodomain and bidomain models suggests that they yield very 
similar results as long as no strong electric fields are applied (1231). To allow comparison with 



(|2ll ). we have used the Courtemanche et al. human atrial model (j24l ) for the ionic currents, with 
the alterations described below. This model is well established and very detailed, taking into 
account all major ionic transport mechanisms and intracellular calcium handling. Note that 
resonant drift in monodomain models was investigated both with atrial (fisi ) and ventricular (0) 
cellular kinetics and results were similar, modulo the difference between steadily rotating and 
meandering spirals. So we expect that our present result should be interesting for ventricular 
fibrillation too, subject to a proper account of other important differences between ventricles 
and atria. 

Numerical methods and parameters 



All simulations were performed by the Cardiac Arrythmia Research Package (CARP) (|26|, |27 



28l ). We used a numerical setup similar to Plank et al. (j2lh . We used a thin sheet of cardiac 
tissue 4 X 2 X 0.02 cm^ with the fibres along the x axis, no-flux boundary conditions and no 
surroundiiig bath. 

As in (j2ll ). the intracellular conductivities in bidomain calculations were fluctuating, with 
conductivities at different points being uncorrelated random numbers within ±50% of the av- 
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Ob j ect / quantity 


Notation / value 


Bidoinain model 


dV 1 


•((T,V(F + 0e)), 






-V • (<T^VV) + le 


iVlUliUU.Ulllcllli lllUClt:! 


C ^- I + 
ot 


• ((T,„Vy) 


Intra- and extra-cellular potentials 


and (pe 


Transmembrane potential 


V = (f)i- 4>e 


Transmembrane ionic current density 


lion, as in (24), with alterations 


External extracellular current density 


Ie{x,y, z,t), by stimulating 


electrodes 


Average cell's surface to volume ratio 


/3 = 1400 cm"^ 


Membrane capacitance per unit area 


Cm = 1.0/iF/cm^ 


Extracellular conductivity, tensor 




. . . , along the fibers 


o'ez = 0.625 S/m 


. . . , across the fibers 


CTey = 0"ez = 0.236 S/m 


Intracellular conductivity, tensor 


<Tj = ai{l + F7]{x,y,z)) 


. . . , fluctuations intensity 


F = 0.5 


. . . , uncorrelated equidistributed noise 


v{x,y,z) G [-1,1] 


. . . , average, along the fibers 


ai^ = 0.174 S/m 


. . . , average, across the fibers 


f^iy = O'iz = 0.019 S/m 


Monodomain tissue conductivity, tensor 




. . . , along the fibres 


amx = 0.146 S/m 


. . . , across the fibres 


CTmy = CTmz = 0.0182 S/m 


Space discretization step 


0.01 cm 


Time discretization steps, most of the time 


10 /is 


. . . , during and 10 ms after single shocks 


1 flS 



Table 1: Details of calculations 



erage, 

aij = aij{l + F7]), (4) 

where j = x, y, z, the level of fluctations was fixed to F = 0.5, and ij G [—1,1] are independent 
equidistributed random numbers. The justification for introducing such fluctuations is that 
although cardiac tissue at a macroscopic scale is frequently approximated as a homogeneous 
bisyncytium, this is not valid at a microscopic scale, and this makes an essential difference 
when an external electric fleld is applied. As the simulation study (j2ll ) showed, the presence of 
fluctuations was fundamental for the mechanism of deflbr illation, although the exact value of 



F was less important. Apart from the F = 0.5 which was the maximal considered in (l2ll ). we 
also tried F = 0.25 in a few test simulations and found the difference not principal. 
Values of the numerous parameters used in this study are presented in Table 1. 

Visualization 



We used CARP's transmembrane voltage V{x, y, t) output to visualize the results of simulations, 
using the colour-coding shown in Fig. [T]C. The same information was used to detect the wave 
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Figure 1: A: Locations of registration electrodes. B: Meandering re-entry with the trajectory 
of the tip (white curve) and the location of the stimulating electrodes. C: Stationary-rotating 
re-entry with the trajectory of the tip. Legend on the right: colour-coding of V{x,y,t). D,E: 
Feedback algorithm with tdeiay = ms (D) and tdeiay = 30 ms (E). Red solid line: registered 
signal, dashed blue line: stimulation shocks, A = Ax 10® //A/cm^. Intervals between shocks are 
the same as the re-entry periods as registered by the electrode. 



tips. The wavefront at time t was defined as the line V{x,y,t) = — 23.75 mV. Then the wave 
tips for that time were defined as the intersections of the front at time t with the front at time 
t — 3 ms. The time delay of t = 3 ms was chosen purely empirically. 

To visualize drift, we used a "stroboscopic" method: we showed positions of the tip syn- 
chronized with the signals detected by the registration electrodes. 

For meandering spirals this required further refinement, since the stroboscopic selection of 
tips produced a congested picture, even without any stimulation: so, a five-petal meandering 
pattern (as shown in Fig. [TJ3) produced five clusters of tips. With stimulation, the tip picture 
becomes even more complicated and unreadable. So we showed only every A^-th stroboscopic 
tip position. Since the meandering patterns were affected by stimulation, we found that the 
optimal value for N was different from 5; we used A = 3. 

A possible alternative to the stroboscopic method is sliding averaging of the trajectories. 
However, it requires careful choice of the averaging window, which may be different in different 
situations, so we found it less convenient. 

Generation of re-entry patterns 

We used two alterations of the Courtemanche human atrial model (0) of the ionic currents, 



(i) To prevent the transmembrane voltages from rising to non-physiological values during 
the defibrillation shocks, an electroporation current was included, see (|2ll ) and refer- 
ences therein. In addition, a formulation for a acetylcholine-dependent potassium current, 
was added (2ll.l29l). This ionic model was used to generate a meandering re-entry. 



I 



K(ACh)^ 

see Fig. [1)3 • 
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(ii) In addition to electroporation and lK{ACh)^ we used a 65% block of the slow inward L-type 
Ca^"*" current coupled with a nine-fold increase in the slow delayed outward K"*" current 
and the rapid delayed outward K"*" current as suggested by Xie et al. (jsol ). This ionic 
model was used to generate a stationary rotating re-entry, see Fig. [TJC. 

To initiate re-entry, we used an S1-S2 protocol. 
Single shock defibrillation benchmark 

Monophasic current shocks Ie{x, y, z, t) were injected into, and withdrawn from, the extracellular 
space via volumes 0.1 x 0.1 x 0.02 cm^, "stimulating electrodes", centered along the left and 
right edges of the slab (Fig. [IB)- The shocks were of rectangular waveform, 5ms duration and 
varied amplitude A. 

Shocks were applied at twelve different timings to, separated by 10ms intervals, after the 
same initial conditions. This covered an entire single rotor cycle, 120ms. A single shock was 
deemed successful if no re-entry was detectable at the moment 500 ms after its end. We define 
the single shock success threshold as the shock amplitude which gives a 50% success rate across 
the twelve timings. 

Resonant drift pacing 

Repetitive low-amplitude shocks of the same waveform and via the same stimulating electrodes 
as above, were applied at the time moments determined by signals received via registration 
electrodes. 

Six different locations for the registration electrodes were used: "point" electrodes 0.02 x 
0.02 X 0.02 cm^ in the top left, top right, bottom left and bottom right corners and "line" 
electrodes (flil ) (of cross-section 0.02 x 0.02 cm^) through the whole medium either horizontally, 
along the fibres or vertically, across the fibres, see Fig. [T]A.. 

The signal from a registration electrode was defined as the average potential at all the nodes 
covered by it. The signal triggered a shock application when it exceeded — 55mV. A shock 
was applied with a delay tdeiay after it had been triggered. In most cases we set tdeiay = 0, see 
Fig. mD). To demonstrate its effect on the direction of the drift, we set tdeiay = 30ms, about a 
quarter of the re-entry period, see Fig. 

In all simulations, the registration electrodes were de-activated for a "blanking time" tbiank = 
50 ms after a shock application. 

For the meandering re-entry, resonant drift pacing was considered successful if the re-entry 
was terminated sooner than the self-termination time of the re-entry, 16000ms. For the sta- 
tionary rotating re-entry, this time was extended to 30000 ms. 

Results 

Single shock defibrillation results 

We have varied the timings, and strength A, of the single shocks to assess the variability of 
the outcomes. Consistent with previous observations (fill ), there were three typical single-shock 
defibrillation outcomes for both the meandering and stationary-rotating patterns: 

• Strong enough shocks annihilated re-entry immediately (Fig. [2]A). 
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Figure 2: Single shock defibrillation, stationary rotating re-entry. A: Termination, A = 16 x 
lO^/uA/cm^ applied at t = 10. B: Breakup, A = U x 10*^ jj. A/ cm^ applied at t = 10. C: 
Displacement, A = 12x 10^ fiA/cui^ applied at t = 10. D: Percentage of successful terminations 
as function of shock amplitude A. Markers: raw data; lines: Bezier approximation. 



• Weaker shocks led to multiple wavebreaks (Fig. [2j3). 

• Weaker still shocks only shift re-entry in space (Fig. [2lC). 

We have found the single shock defibrillation threshold to be ^ = 14 x 10^ /iA/cm^ for the 
meandering patterns and 18 x 10^/iA/cm^ for the sationary-rotating patterns (Fig. [2p). 

Resonant drift pacing results 

We tested resonant drift pacing on the meandering, and stationary rotating, re-entrant pat- 
terns with shock amplitudes lower than the corresponding single shock defibrillation success 
thresholds. 



Meandering re-entry 

Point registration electrodes. Depending on the shock strength, the following outcomes 
were observed. 
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t = 111 ms, n = 1 t = 2253 ms, n = 1 




i = 3715 ms, n = 30 t = 4087 ms,n = 33 





t = 4835 ms, n = 42 t = 5000 ms, n = 43 




i = 118 ms, n = 1 t = 2252 ms, n = 18 t = 2990 ms, n = 24 




t = 3834 ms, n = 30 t = 4982 ms, n = 39 t = 5625 ms, n = 45 




t = 6247 ms, n = 51 i = 6419 ms, n = 52 i = 6460 ms, n = 53 



A B 

Figure 3: Resonant drift pacing of a meandering re-entry. A: A = 1 x 10^ /xA/cm"^, re-entry 
drifts to the boundary and terminates. Here and below, cross: location of the registration 
electrode, "n = . . .": number of shocks applied so far, white dots: positions of the spiral tip 
at every third registered period. B: A = 2 x 10^/uA/cm^, original re-entry "A" drifts to the 
boundary and terminates. Additional re-entrant patterns "B", "C" and "D" are generated by 
the shocks and trigger further shocks themselves. "B" terminates at the boundary, "C" and 
"D" annihilate each other. NB: resonant drift pacing handles multiple patterns. 



• The original re-entry drifted, roughly towards the registration electrode, until termination 
on an inexcitable boundary, see Fig. [3]A. 

• Secondary re-entrant patterns were generated near stimulating electrodes; both the origi- 
nal and newly generated re-entries drifted. An example is shown in Fig. [3l3. The original 
re-entry "A" was terminated by the first 24 shocks. However, a secondary re-entrant 
pattern "C" was generated by a shock in the process. After "A" terminated, it was "C" 
that triggered the subsequent shocks. A third re-entrant pattern, "B" was generated by a 
further shock. Due to the proximity of "B" to the registration electrode, it took over the 
control of the shock applications. Re-entry "B" drifted to the boundary and terminated 
after a further 15 shocks. Then "C" was again the solitary re-entry and took over the 
control. A fourth re-entry, "D" was generated by a shock and, due to its proximity to the 
registration electrode, took over the control. After 3 shocks, "D" collided with "C" and 
they annihilated each other without reaching an inexcitable boundary. 

Note that mutual annihilation of re-entrant patterns in the bulk of the tissue could be a 
plausible explanation why Pak et al. (jsil ) saw no evidence of annihilation at the ventricular 
boundaries when testing a protocol similar to the resonant drift pacing. 



Line registration electrodes 

• Horizontal electrode: With the registration electrode along the fibres, the re-entry drifted 
at a small angle to it (Fig. |3]A) . 

• Vertical electrode: With the registration electrode across the fibres, the drift was usually 
more complicated. An example is shown in Fig. |3J3. The re-entry initially drifts upwards 
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Figure 4: Resonant drift pacing with line registering electrodes. A: Horizontal electrode 
along the fibres, A = 2 x 10^/iA/cm^. The re-entry drifts resonantly at an angle to the 
electrode, reaches the boundary and terminates. B: Vertical electrode across the fibres, 
A = 1 X 10^ /iA/cm"^. The re-entry drifts upwards for 21 shocks, then changes direction and 
reaches the boundary after a further 109 shocks. 



along the electrode. After 21 shocks, the drift turns downwards, until it reaches an 
inexcitable boundary after a further 109 shocks and terminates. The trajectory of the 
drift crosses itself, which is not allowed in the asymptotic theory of drift of rigidly rotating 



spiral waves (llTI ) and was not observed in our simulations of a stationary rotating re-entry. 



so we consider it a new feature due to meander. 

The time taken for termination of all re-entrant activity using different locations of the 
registration electrode is shown in Fig. OA. It shows that resonant drift pacing can successfully 
terminate meandering re-entrant patterns using shocks 14 times weaker than the single shock 
defibrillation threshold. Termination was achieved in all simulations in this series. 

Among the registration electrodes, the vertical line was the most successful, as it produced 
the largest proportion of successful terminations in the fastest time except at the smallest 
amplitudes. Using a point-electrode in the top right location was the least successful. 

We conclude that the termination time depends on the mutual position of the registration 
electrode, anode, cathode, initial position of the re-entry, tissue size and fibre orientation. 

For amplitude ^4 = 1 x 10^/iA/cm^, re-entry termination time for the meandering pattern, 
in most cases, was longer than its self-termination time 16000ms. Thus, we did not consider 
lower values for the stimulation amplitudes. 



Stationary rotating re-entry 

Fig. shows the times taken for annihilation of all re-entrant waves for the different locations 
of the registration electrode. Resonant drift pacing can terminate stationary-rotating patterns 
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Figure 5: Resonant drift pacing results. A: Meandering re-entry. Left: termination time as a 
function of shock amplitude smooth curve: Bezier approximation. Right: the raw data, time 
in ms and amplitude A in lO^/uA/cm'^. B: Same, for a stationary rotating re-entry, oo: infinite 
loop. 



using amplitudes 18 times lower than the single shock threshold. The combination of shock 
strength with the location of the registration electrode affects the probability of success and the 
time taken for termination. 

Termination of stationary rotating re-entry was not always observed within the 30000 ms 
limit. Sometimes the algorithm gets caught in an infinite loop (Fig. [6|), when the original re- 
entrant pattern has been terminated, but shocks produce new wavefronts which trigger further 
stimulation producing further wavefronts and so on. 



Direction of drift 



The theory of the resonant drift (jl5l ) predicts that the drift direction depends on the stimulation 
phase. This dependence leads to the relationship between the drift direction and the delay in 
the case of a feedback controlling forcing (17], |32| ). A delay by a certain fraction of a spiral's 
period causes a change in the drift direction by the same fraction of 360° . We have verified that 
it works in this realistically anisotropic bidomain model as well. 

Fig. [7]A shows fragments of two such simulations, with the same initial conditions and 
different values of tdeiay, and 30 ms, which is about a quarter of the re-entry period. The 
direction of the drift in these two cases differs roughly by about 90°, with account of the 
anisotropy, in agreement with (17, [13). 
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Figure 6: Infinite loops in resonant drift pacing. A: ^ = 10 x 10^/iA/cm^. New waves with 
breaks are initiated by shocks and trigger further shocks. B: ^ = 2 x 10^/xA/cm^. New waves 
without breaks are initiated by shocks and trigger further shocks. 



Using a time delay to induce success 

We have mentioned in section that resonant drift pacing may fail through an infinite loop of 
wavefront eliminations and creations. If newly created wavefronts are unbroken as in Fig. [6(3 
then the re-entrant activity will vanish if the stimulation is stopped at any time, i.e. the success 
is achieved. 

A genuine failure happens when newly created waves are broken and form new re-entries, as 
in Fig. [B]A. The chances of getting in such a loop depend on the drift trajectory, which depends 
on the location of the registration and stimulation electrodes and the feedback delay. Hence such 
loops should be avoidable by an appropriate choice of the electrode locations and/or of tdeiay 
In a real-life situation, such choice would require plenty of information and may be difficult to 
make. Not all theoretical solutions {e.g. optimal positions of the electrodes) may be possible to 
implement in practice. Further, the optimal parameters may vary from one arrhythmia episode 
to another so difficult to predict a priori. 

Of all the mentioned parameters, tdeiay is the easiest to change, and it may even be adjusted 
in real time, during pacing. Hence we propose that a change of the feedback delay may be used 
to discontinue an infinite loop even after it has already started. 

Fig. [7p illustrates the feasibility of this approach. It is a re-run of the simulation of Fig. [6]A., 
with tdeiay increased to 30 ms at t = 1000ms. This broke the loop and terminated the re-entry. 
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i = 976 ms, n = 8 
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i = 280 ms, n = 3 t= 540 ms, n = 7 




i = 1000 ms, n = 14 i = 1222 ms, n = 16 




t = 1304 ms, n = 17 i = 1476 ms, n = 1 




t = 1854 ms, n = 20 t = 1936 ms, n = 21 




i = 2114 ms, n = 22 t = 2490 ms, n = 24 




^ t = 2572 ms, n = 25 t = 2800 ms, n = 26 

Figure 7: Resonant drift of a meandering pattern with t^eZaj/ = ms (A) and tdeZaj/ = 30 ms, 
^ = 1 X lO^^A/cm'^ (B). C: Tip trajectories of A, B enlarged. D: Effect of increasing tdeiay to 
30 ms aXt = 1000 ms to the infinite loop shown in Fig. [6]A: the re-entry is extinguished. Arrows: 
generation of new waves and the direction of drift. NB: resonant drift pacing handles multiple 
patterns 



Discussion 

Summary of results: Resonant drift pacing works in a realistic model of elec- 
tric field action on cardiac tissue 

This simulation study presents results which can be used as experimentally testable hypotheses. 
It has been based on the understanding of the mechanism by which external electric current 
affects excitation and propagation via non-uniformity of the electric field and the heterogeneity 
of the tissue conductivities. This mechanism implies both a strong action near electrodes due 
to the electric field inhomogeneity, and a weaker, but far reaching action in the bulk of the 
tissue due to omnipresent tissue heterogeneities. This model has been used before to simulate 
single-shock defibrillation (fill ). Here we have demonstrated that the far- field action of the 
electric current can cause, under repetitive stimulation, a drift and elimination of the re-entrant 
sources. In our simulations, resonant drift pacing can eliminate re-entry, including multiple 
re-entries, with high probability and within acceptable time, at amplitudes much lower than 
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single shock defibrillation. If we allow 10 s for low- voltage re-entry termination, as is the case 
for antitachycardia pacing, the required shock strength with resonant drift pacing is 12 to 15 
times smaller than with a single shock. 

Our simulations have shown that a major possible obstacle to elimination of re-entrant 
sources is the possibility that the electric shocks create new re-entrant waves while eliminating 
existing ones, which may lead to infinite loops of annihilation and creation. The creation of the 
new sources occurs near the shock electrodes where the electric field is highly inhomogeneous, 
and the 12 to 15-fold decrease in defibrillation threshold is observed notwithstanding this effect. 
The highly inhomogeneous electric field is created by point electrodes, so bigger electrodes, 
which create more homogeneous fields, should perform even better. 

Our simulations have also demonstrated yet another way to overcome infinite annihilation- 
creation loops. This is due to the dependence of the direction of resonant drift on the time 
delay in the feedback loop. This dependence has been predicted and observed in isotropic 
monodomain models with simplified description of the electric field action 
Here we demonstrated that this dependence is still observed in the present more realistic model 
and, moreover, it can be used to discontinue an infinite annihilation-creation loop, via a change 
in the feedback delay once such a loop has been detected. 

Although we used atrial tissue kinetics, there are indications that the exact sort of excitable 
kinetics is not too important for the properties of resonant drift (jl9l . l25l). So we expect that our 
present results may be interesting for ventricular fibrillation, too. 



Comparison with previous research 

Anti-tachycardia pacing In clinical practice, antitachycardia pacing is as efficient and as 
safe as single shock defibrillation, even for fast ventricular tachycardias with up to 250 beats/min (fl3 ). 
The exact mechanisms underlying success or failure of antitachycardia pacing are widely unex- 
plored. The traditional understanding implies stimulation which is faster than the anatomical 
re-entry and which engages all the larger area of tissue until it reaches an "isthmus" of the re- 
entry and blocks it. It is conceivable, however, that in some cases antitachycardia pacing may 
be unwittingly applied to a functional rather than anatomical re-entry and then its mechanism 
could be the resonant drift. In such cases, feedback-controlled resonant drift pacing will do a 
better job than pacing with prescribed frequency used in antitachycardia pacing. Although the 
probability of antitachycardia pacing success is high, when it fails a single strong shock has to 
be applied. 



Unpinning It has been speculated that anatomically or functionally anchored ventricular 
tachycardias are less likely to be terminated by antitachycardia pacing. For these cases, it has 
been suggested that a weak shock, with an amplitude in the range used in this study, can be 
applied to unpin the re-entry from the obstacle by a virtual electrode polarization mechanism 



34l ). Such a stimulus must be correctly timed, which can be achieved by synchronizing 



it to the signal from a registration electrode, with a correctly chosen delay with respect to 
that signal, which in practice may require scanning through possible delays. The resonant drift 
pacing also requires synchronization with a registered signal, and we have seen that change in 
the delay during the course of pacing may also be beneficial. This makes the unpinning protocol 
operationally very similar to our resonant drift pacing protocol, and in any particular case it 
may not be possible to say with certainty which mechanism has worked. 
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Experiments with feedback-controlled pacing Pak et al. (j3ll ) eliminated ventricular fib- 
rillation in rabbit by multisite pacing synchronized with optical signal from a fixed reference 
site. Their stimulation protocol was similar to resonant drift pacing considered here, except 
they used more pacing electrodes and aimed to deliver shocks only when pacing site was in 
an excitable gap, under the implicit assumption that shock delivered to a site in an absolute 
refractory state could not possibly affect the re-entry, i.e. a priori ignoring far-field effects. 
Despite this artificial self-limitation, they were successful at defibrillating with shock strengths 
an order of magnitude lower than single shock defibrillation and with better success rate than 
overdrive and high-frequency pacing. It was argued in (jsil ) that earlier mathematical modeling 



of resonant drift only demonstrated termination of a single re-entrant wave (see, however, ()32l )) 
and their experiments did not show evidence of extinction of re-entries on ventricular borders. 
Our present simulations indicate that 

• far-field effects can make it worthwhile issuing shocks regardless of the tissue state at the 
pacing site, 

• the resonant drift pacing can work for multiple re-entrant waves (providing that the dif- 
ficulties of new wavefronts being initiated by the shocks are overcome), and 



annihilation at the boundaries is not a necessary feature of resonant drift elimination, as 
re-entry sources can annihilate with each other. 



Thus, we find the results by Pak et al, contrary to their own interpretation, as strong experimen- 
tal evidence of the resonant drift pacing being an efficient method of low voltage defibrillation. 



Clinical implications 

Resonant drift pacing may present an alternative or a supplement for existing (single shock 
defibrillation, synchronized cardioversion, antitachycardia pacing, overdrive pacing, and high- 
frequency pacing) and proposed (unpinning) therapies. The theoretical mechanism underlying 
this protocol is for functionally determined re-entries, i.e. probably for higher frequency tachy- 
cardias and fibrillations. Due to similarity between resonant drift and unpinning pacing proto- 
cols, the two protocols may be possible to combine into one, which would work via unpinning 
for pinned (monomorphic tachycardias), and via resonant drift for unpinned re-entrant waves 
(polymorphic tachycardias and fibrillations). 

There is evidence that some atrial arrhythmias are not re-entrant, but are due to rapid 
ectopic focal activity (jssl ). The proposed method might work on such arrhythmias via a com- 
pletely different mechanism, say overdrive suppression, but this can only be clear after further 
studies. 

As the study presented here is purely theoretical, the key question is whether the low-energy 
approach to defibrillation considered here could work in real heart. Clearly this question can 
only be answered by experimental studies, and there is clear scope to advance the field by 
interaction between modeling studies such as this one and experimental studies. 



Study limitations 

Apart from inevitable limitations inherent in mathematical modeling compared to experimental 
study, which are due to limited current knowledge, we have made a number of simplifications. 
Our model tissue was essentially two- rather than three-dimensional. Its geometry was not 
realistic (rectangular shape). The model tissue lacked any macroscopic inhomogeneities, such 
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as trasmural, center-periphery or base-apex gradient of excitability properties, or variations of 
the conductivity tensor. The position and geometry of the electrodes was arbitrary. Values 
of some of the parameters lack reliable experimental foundation, including some of the most 
important ones, such as F, the intracellular conductivity fluctuations amplitude. The spatially 
uncorrelated structure of conductivity fluctuations is an idealization. These limitations can be 
lifted as more experimental data become available and via further, more detailed, simulation 
studies. 

Resonant drift pacing may potentially work for ventricular tachycardia and fibrillation too. 
However, the ventricles have thick walls and overall more complicated anatomy, which may 
pose extra difficulties for defibrillation. So any extrapolation of our present results to ventricles 
should be thought of with caution, bearing in mind possible complications. There are modeling 
studies suggesting that at least some of the new difficulties arising in resonant pacing of three- 
dimensional re-entries can be successfully overcome (see (0) and references therein), however 
those were in overly simplified models of exciable media, and more realistic computational 
studies are still desirable. 

Conclusion 

Our results show that in a bidomain model of cardiac tissue, with microscopic conductivity 
fluctuations, re-entrant waves can be annihilated using feedback-controlled repetitive stimula- 
tions by inducing resonant drift and directing the re-entry towards an inexcitable boundary. 
If associated difficulties are overcome, termination using this approach is achieved, with high 
probability and within acceptable time, at a fraction of the conventional shock strength. The 
direction of the drift can be managed by choosing the location of the electrodes and the time 
delay of the shock application. 

Difficulties occur due to new wavefronts being initiated from the shock electrodes. However, 
numerical simulations allow a detailed insight into this problem and can be used as a tool to 
suggest solutions. 

Our findings present experimentally testable hypotheses for what we expect to observe in 
real cardiac tissue. There is scope to advance the field of low-energy defibrillation by interaction 
between modeling studies such as this and experimental studies. 
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